Carrier Density Collapse and Colossal Magnetoresistance in Doped Manganites 
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A novel ferromagnetic transition, accompanied by carrier density collapse, is found in doped charge- 
transfer insulators with strong electron-phonon coupling. The transition is driven by an exchange 
interaction of polaronic carriers with localized spins; the strength of the interaction determines 
whether the transition is first or second order. A giant drop in the number of current carriers during 
the transition, which is a consequence of bound pairs formation in the paramagnetic phase close 
to the transition, is extremely sensitive to an external magnetic field. This carrier density collapse 
describes the resistivity peak and the colossal magnetoresistance of doped manganites. 
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The interplay of the electron-phonon and exchange in- 
teractions (l]-|j| is thought to be responsible for many 
exotic properties of oxides ranging from high-T c super- 
conductivity in cuprates p| to colossal magnetoresis- 
tance (CMR) and ferromagnetism in doped manganites 
|^|-^|. A huge negative magnetoresistance was observed 
in doped perovskite manganites Lai-^D^MnOa (D=Ca, 
Sr, Ba) close to the ferromagnetic transition in a certain 
range of doping x w 0.15 — 0.4 0-0], and this raised a 
question of possible applications. 

The metal-insulator transition in lanthanum mangan- 
ites has long been thought to be the consequence of a 
double exchange (DEX), which results in a varying band 
width of holes doped into the Mn 3+ d-shell |10|, as func- 
tion of the doping concentration and temperature. Re- 
cently it has been realized pj] |, however, that the effec- 
tive carrier-spin interaction in DEX model is too weak to 
lead to a significant reduction of the electron band width 
and, therefore, cannot account for the observed scatter- 
ing rate fl| (see also Ref. ||) or for localization induced 
by slowly fluctuating spin configurations flli| . In view of 
this problem, it has been suggested JllJ that the essential 
physics of perovskite manganites lies in the strong cou- 
pling of carriers to Jahn- Teller lattice distortions. The 
argument Jll| was that in the high-temperature state 
the electron-phonon coupling constant A is large (so that 
the carriers are polarons) while the growing ferromag- 
netic order increases the bandwidth and thus decreases 
A sufficiently for metallic behavior to set in below the 
Curie temperature T c . A giant isotope effect ||l5| , the 
sign anomaly of the Hall effect, and the Arrhenius be- 
havior of the drift and Hall mobilities jl6| over a temper- 
ature range from 2T C to 4T C unambiguously confirmed 
the polaronic nature of carriers in manganites. 

However, an early established unusual relation be- 
tween magnetization and transport below T c have led 
to a conclusion that the polaronic hopping is the preva- 
lent conduction mechanism also below T r [|17L Low- 



temperature optical ||l8|-p0||, electron-energy- loss (EELS) 
ifjlf and photoemission spectroscopies [^2| showed that 
the idea of a 'metalization' of manganites below 

T c is not tenable. A broad incoherent spectral feature 
|18 -p0tp2[ and a pseudogap in the excitation spectrum 
[22,g3[ were observed while the coherent Drude weight 
appeared to be two orders of magnitude smaller jl9| than 
is expected for a metal. EELS |2l| confirmed that man- 
ganites are charge-transfer type doped insulators having 
p-holes as the current carriers rather than d (Mn 3+ ) elec- 
trons. The photoemission and O Is x-ray absorption 
spectroscopy of Lai-^Sr^MnOs showed that the itiner- 
ant holes doped into LaMnC>3 are indeed of oxygen p 
character, and their coupling with d local moments on 
Mn 3+ ions aligns the moments ferromagnetically 
Moreover, measurements of the mobility |^,^5| do not 
show any field dependence. The calculated resistivity is 
in poor agreement with the data and the characteristic 
theoretical field (~ 15T) for CMR is too high compared 
with the experimental one (~ 4T) |0]. As a result, self- 
trapping above T c and the idea of metalization below T c 
do not explain CMR either. Carriers retain their pola- 
ronic character well below T c , as manifested also in the 
measurements of resistivity and thermoelectric power un- 
der pressure 

In the present paper, we propose a theory of the 
ferromagnetic/paramagnetic phase transition in doped 
charge-transfer magnetic insulators accompanied by a 
current carrier density collapse (CCDC) and CMR. Tak- 
ing into account a tendency of polarons to form bound 
pairs and the (competing with binding) exchange inter- 
action of p polaronic holes with d electrons, we find a 
novel ferromagnetic transition driven by non-degenerate 
polarons. As a result we describe the magnetization 
and temperature/field dependence of the resistivity of 
Lai-^CaaMnOs close to T c in a region 0.15 < x < 0.4. 

The Hamiltonian containing the physics compatible 
with the experimental observations is 
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where E^ is the LDA energy dispersion [g7| j ^ks is the an- 
nihilation hole operator of a (degenerate) p oxygen band 
with spins s =f and |, J p( 2 is the exchange interaction 
of p holes with four d electrons of the Mn 3+ ion at the 
site i, mk = h^hy^ — /iL/ikj., S* is the z-component 
of Mn 3+ spin, which is 5 = 2 due to the strong Hund 
coupling, WhuikIj °f the four d-electrons on Mn 3+ sites, 
N is the number of unit cells. The two last terms of 
the Hamiltonian describe the coupling of p holes with 
phonons and the phonon energy respectively (w q is the 
phonon frequency, 7 q = ^ 7- q is the coupling constant 
If the holes were doped into d shell instead of p 
shell, the Hamiltonian would be similar to the Holstein 
t-J model H with about the same physics of CMR as 
proposed below. 

Essential results are readily obtained within the 
Hartree-Fock approach for the exchange interaction |2S|] 
and the Lang-Firsov polaron transformation, which elim- 
inates terms linear in the electron-phonon interaction in 
Eq- (0) Pi- Thus, we find spin-polarized p bands 
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where £k = E^e^ 9 , and g 2 ~ j 2 describes the pola- 
ronic band narrowing and the isotope effect fl5|, a is 
the normalized thermal average of the Mn spin (3(); H 
is the external magnetic field, and \ib is the Bohr mag- 
neton. The p — d exchange interaction depends only on 
total (average) magnetization because we assume that 
the system is homogeneous. The magnetization of Mn 3+ 
ions is given by 

— /o*\/o d ( J P dm + 2g Mn ^ B H\ 
*=(S t )/S = B s ^ _ j, (3) 

with m the absolute value of the magnetization of holes 
determined as 



m = 
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Here B s (x) = [1 + 1/(25)] coth[(5 + l/2)x] - 
[1/(25)] coth(a;/2) is the Brillouin function, gMn the 
Lande g- factor for Mn 3+ in a manganite, N p (e) the den- 
sity of states in the narrow polaron band, and f p (e^ s ) — 
[y _1 exp(eks/fcsr) + 1] _1 the Fermi-Dirac distribution 
function with y = exp(/i/fcsT) determined by the chemi- 
cal potential \x. Note that for J p d < (antiferromagnetic 



coupling) the main system of equations remains 
the same after a substitution J p d — > \J P d\- 

Along with the band narrowing effect, the strong 
electron-phonon interaction binds two holes into a pair 
(bipolaron) |jj . The bipolarons are practically immobile 
in cubic manganites because the electron-phonon interac- 
tion is too strong in contrast with cuprates, where bipo- 
larons are mobile owing to their geometry and a moderate 
coupling with phonons | p9[ . 

If these bound pairs are extremely local objects, two 
holes on the same oxygen, then they will form singlet. 
If, however, these holes are localized on different oxy- 
gens, then they may form a triplet state. Because of 
their zero spin, the only role of singlet bipolarons in 
manganites is to determine the chemical potential /z, 
which can be found with the use of the total carrier 
density per cell x as 2 J dENb p (E)fb p (E) = x — n, 
where Nh p (E) the density of bipolaronic states, fb P (E) = 
{y~ 2 exp[(i? — A)/fcsT] — 1} _1 the bipolaron distribution 
function, A the bipolaron binding energy, n is the den- 
sity of single (unbound) hole polarons, which are the only 
current carriers in manganites, and x is the doping con- 
centration. 

It is the localization of p-holes into immobile bound 
pairs combined with their exchange interaction with the 
Mn d A local moments that are responsible for CMR. The 
density of these pairs has a sharp peak at a ferromag- 
netic transition when system is cooled down through the 
critical temperature T c . Below T c the binding of po- 
larons into pairs competes with the ferromagnetic ex- 
change which tends to align the polaron moments and, 
therefore, breaks those pairs apart. These competing in- 
teractions lead to unusual behavior of CMR materials 
and a huge sensitivity of their transport to external field. 

To illustrate the point we assume that T c is compa- 
rable with the polaron, W, and bipolaron band widths 
IpOj . Then (bi)polarons are not degenerate in the rel- 
evant temperature range, f p ~ yexp(— E/ksT) and 
/bp ~ y 2 exp[(A — EJ/kgT], and we can evaluate in- 
tegrals reducing the system of mean field equations to 
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Here we use the dimensionless temperature t = 
2ksT '/ '(JpdS), magnetic field h — 2[LbH j '(J p dS), and the 
binding energy S = A/(J p( j5), while v{= 3) is the degen- 
eracy of the p band. 

The polaron density n is determined by Eq. (|B|) with 
(7 = above T c . At the critical temperature, the polaron 
density has a minimal value n c ~ (2a;) 1 / 2 exp(— 8/t c ), it 
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then grows exponentially with temperature and saturates 
at n = (1 + 2.T) 1 / 2 — 1. The remarkable observation is that 
there is a sharp increase of the polaron density at tem- 
peratures below T c . The physical origin of the unusual 
minimum of the current carrier density at T c lies in the 
instability of bipolarons below T c due to the exchange in- 
teraction of polarons with d electrons. The spin-polarized 
polaron band falls below the bipolaron band with de- 
creasing temperature, so that all carriers are unpaired at 
T = if J pd S > A. 

Linearizing (^j)-(f^) we find the critical temperature in 
zero magnetic field is t c = (nc/2) 1 / 2 , where the polaron 
density at the transition n c is determined by 

nl' 2 \n 2{ ~ X ~ nc) =2*' 2 5. (9) 
n 2 . 

This equation has solutions only for S below some critical 
value S c (x), Fig. l(b, inset). The numerical solution of 
the system Eqs. (||)-(^) shows that for S > 8 c (x) the 
ferromagnetic phase transition is first order with jumps of 
the polaron density and the magnetization, as observed 
plfl , Fig. 1(a). The transition is continuous when 5 < 
S c (x), Fig. 1(b). A relatively weak magnetic field has a 
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FIG. 2. Inverse polaron density for different magnetic 
fields, A/ JpdS = 0.5, doping x — 0.25. Note that the tran- 
sition is a strong first order, and becomes continuous only 
when the external magnetic field exceeds some critical value. 
Inset: schematic of polaron (P) and bipolaron (BP) densities 
of states at temperatures below and above T c for up (|) and 
down (J.) spin moments. The pairs (BP) break below T c if 
exchange J p dS exceeds the pair binding energy A. 




FIG. 1. Polaron density n/x and magnetization 

as a function of temperature t — 2ksT/( J p dS) in a 
doped charge-transfer insulator near (a) the first order, 
8 = A/(JpdS) = 0.5, and (b) second order, 5 = 0.25, phase 
transitions (doping x — 0.25). Inset: the critical value of 
the relative binding energy of polaron pairs 5 separating the 
regions of the first and second order phase transitions. 



drastic effect on the inverse carrier density, 1/n, near the 
first order transition, as shown in Fig. 2. As a result, the 
resistivity p = l/(enfi p ) has a sharp maximum, which is 
extremely sensitive to the magnetic field in the vicinity 
of T c . 

It is assumed, as is usually the case, that the triplet 
states always lie higher in energy than the singlet state. 
If the singlet-triplet separation becomes smaller than the 
gap, J st < A, then, because of a higher number of 
the triplet states, their thermal population leads to a 
deeper minimum in the density of polarons. We make 
an essential assumption that the exchange between spins 
on Mn and triplet bipolarons is suppressed because the 
bipolarons are strongly localized. Otherwise, the triplet 
bound pairs, if they were formed in the paramagnetic 
phase, can survive in the ferromagnetic phase thus re- 
ducing the carrier density collapse. One can draw an 
analogy of this situation with singlet magnetism, e.g. in 
Pr compounds p^ ]. 

In fact, our theory, Eqs. (^)-(^), describes all the major 
features of the temperature/field dependence of p(T) [Q, 
with a temperature and field independent polaron drift 
mobility fi p p3[ | in the experimental range of the mag- 
netic field, Fig. 3. That suggests that CCDC is the 
origin of CMR. In general, one has to take into account 
the temperature dependence of the polaron mobility to 
extend our theory for temperatures far away from the 
transition. 

We have also compared this scenario with the localiza- 
tion of p- holes due to a random field with a gap A/2 be- 
tween localized impurity levels and the conduction band 
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FIG. 3. Resistivity of Lao.7sCao.2sMn03 calculated 
within the present theory for the pair binding energy 
A = 900K, and polaron exchange with the localized Mn 3+ 
spins JpdS — 2250K (a), compared with experiment [7] on 
panel (b) . Note an extreme sensitivity of the theoretical resis- 
tivity to external magnetic field (a), observed experimentally 
(b) for the doped manganite. 



p4| . We have found the same features of the phase tran- 
sition in zero field. However, Zeeman splitting of the 
impurity states, and a different behavior of y with den- 
sity, makes the transition far less sensitive to the mag- 
netic field. As a result, no quantitative description of the 
experimental CMR has been found in this case. 

In conclusion, we have found that a few non-degenerate 
polarons in the p band polarize localized d electrons be- 
cause of a huge density of states in the narrow polaronic 
band. For a sufficiently large p — d exchange J pc iS > A, 
we have obtained current carrier density collapse at the 
transition owing to the formation of immobile bipolarons 
in the paramagnetic phase with the binding energy A 
x [[p. Competition between the binding energy of po- 
larons, which promotes a formation of local pairs, and 
their exchange interaction with d electrons, which breaks 
them at lower temperatures, results in a huge negative 
magnetoresistance close to the ferromagnetic transition. 

We have explained the resistivity peak and the colos- 
sal magnetoresistance of doped perovskite manganites, 
Fig. 3, as a result of the current carrier density collapse. 
Depending on the ratio A/(J p dS), the transition is first 
or second order, Fig. 1. It is not clear at present whether 
the main idea underlying this picture, the assumption of 
singlet bound states of charge carriers, is true for man- 
ganites. The available experimental data, e.g. the tun- 
neling spectroscopy |2^| suggests that it is. Our goal is to 
stimulate a wider discussion and new experiments in this 
direction. We expect that the present theory is general 
enough to also account for the giant magnetoresistance 



observed in pyrochlore manganites M . 

We acknowledge useful discussions with A.R. Bishop, 
G.A.D. Briggs, D.M. Edwards, M.F. Hundley, P.B. Lit- 
tlewood, S. von Molnar, V.G. Orlov, W.E. Pickett, 
D.J. Singh, S.A. Trugman, and R.S. Williams. We thank 
G. Aeppli, D.S. Dessau, H.T. Kim, A. P. Ramirez, and 
G.-m. Zhao for the information on their data and help- 
ful discussions. 



[3] 

[5] 

[6] 
[7] 
[8] 
[9] 
[10] 



[11] 

[12] 
[13] 

[14] 
[15] 
[16] 
[17] 

[18] 
[19] 
[20] 
[21] 
[22] 
[23] 

[24] 
[25] 
[26] 
[27] 

[28] 

[29] 



Email asa21@cus.cam.ac.uk 
Email alexb@hpl.hp.com 

A.S. Alexandrov and N.F. Mott, Polarons and Bipolarons 
(World Scientific, Singapore, 1995). 

A.R. Bishop and M.I. Salkola, in Polarons and Bipo- 
larons in High-T c Superconductors and Related Materi- 
als, edited by E.K.H. Salje et al. (Cambridge Univ. Press, 
Cambridge, 1995), p. 353. 

H. Fehske et al, Phys. Rev. B 51, 16582 (1995). 
J.G. Bednorz and K.A. Miiller, Z. Phys. B64, 189 (1986). 
J.H. Van Santen and G.H. Jonker, Physica 16, 599 
(1950). 

S. Jin et al, Science 264, 413 (1994). 
P. Schiffer et al, Phys. Rev. Lett. 75, 3336 (1995). 
Y. Shimakawa et al, Nature (London) 379, 53 (1996). 
A. P. Ramirez, J. Phys. Condens. Matter 9, 8171 (1997). 

C. Zener, Phys. Rev. 82, 403 (1951); P.W. Anderson 
and H. Hasegawa, Phys. Rev. 100, 675 (1955); P.G. dc 
Gennes, Phys. Rev. 118, 141 (1960). 

A.J. Millis et al, Phys. Rev. Lett. 74, 5144 (1995); A.J. 
Millis et al, ibid. 77, 175 (1996). 

D. M. Edwards (to be published). 

M.E. Fisher and J.S. Langer, Phys. Rev. Lett. 20, 665 
(1968). 

CM. Varma, Phys. Rev. B54, 7328 (1996). 

G. -m. Zhao et al, Nature (London) 381, 676 (1996). 
M. Jaime et al, Phys. Rev. Lett. 78, 951 (1997). 

M.F. Hundley et al. Appl. Phys. Lett. 67, 860 (1995); 
R.H. Heffner et al. Phys. Rev. Lett. 77, 1869 (1996). 
Y. Okimoto et al, Phys. Rev. Lett. 75, 109 (1995). 
Y. Okimoto et al, Phys. Rev. B 55, 4206 (1997). 
K.H. Kim et al, Phys. Rev. Lett. 81, 1517 (1998). 

H. L. Ju et al, Phys. Rev. Lett. 79, 3230 (1997). 
D.S. Dessau et al. Phys. Rev. Lett. 81, 192 (1998). 
A. Biswas et al, bond-mat/9806084 J.Y.T. Wei et al, 



Phys. Rev. Lett. 79, 5150 (1997) 
T. Saitoh et al, Phys. Rev. B 51, 13942 (1995). 
H.-T. Kim et al. (to be published). 
J.-S. Zhou et al, Phys. Rev. Lett. 79, 3234 (1997). 
W.E. Pickett and D.J. Singh, Phys. Rev. B 53, 1146 
(1996). 

See, e.g.: K. Yosida, Theory of Magnetism (Springer- 
Verlag, Berlin, 1996), p. 201. 

Objections, raised by V.K. Chakraverty et al. [Phys. 
Rev. Lett. 81, 433 (1998)] w.r.t. the bipolaron model 
of cuprates, were shown to be completely erroneous (see 



4 



A.S. Alexandrov, |cond-mat/9807185| ). 
[30] Otherwise there will be an inessential factor (T/W) d ^ 2 in 

Eq. (^) for n, where d is the dimensionality of a system. 
[31] H. Kuwahara et al, Science 270, 961 (1995). 
[32] J. Mrachkov and E. Leyarovski, Z. Phys. B 82, 233 

(1991). 

[33] With our low polaron density at the transition the po- 
laron mobility is fi p = 0.2 cm 2 /Vs for x = 0.33 Jl6| , 
and about 0.03 cm 2 /Vs for x = 0.25 that lies in the 
range typical of polaronic conductors like Ti02 at room 
temperature 0. 

[34] A.S. Alexandrov and A.M. Bratkovsky (to be published). 



5 



